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I ABSTRACT: We examine the dependence of decoherence on the spec- 
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tral density of the environment as well as on the initial state of the system. 
We use two simple examples to illustrate some important effects. 
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O ' Introduction 

o ■ 

O ' Decoherence plays a major role in the transition from quantum to clas- 
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sical and has attracted much attention in recent years (see Zurek (1991)). 
The analysis of this process may allow us to understand in detail the mecha- 
nism that prevents observations of some quantum systems in superpositions 



^ ! of macroscopically distinguishable states. In the light of new technologies 

it can also help us to devise experiments to probe the fuzzy boundary be- 
tween the quantum and the classical world. The interaction with an exter- 
nal environment is the mechanism responsible for the supression of quantum 
interference effects. Therefore, there are several questions that arise natu- 
rally: How dependent on the environment decoherence is? What are the 
time scales involved in this process? How are some preferred states of the 
system dynamically chosen? In this paper we will report on recent work 
where some of these questions are addressed. 



As a first point, let us clarify what we mean here by decoherence. Within 
the Gell-Mann and Hartle version of the consistent histories formulation of 
quantum mechanics, based on earlier work by Griffiths and Omnes, (see 
contributions in this proceedings) the term decoherence is used in place of 
"consistency" which is the condition that, if satisfied, allows us to assign 
probabilities to members of sets of coarse grained histories of a closed sys- 
tem. The Decoherence Functional is the basic diagnostic tool used in this 
framework. On the other hand, in previous works originated in quantum 
measurement theory, a different notion of decoherence was used. Measure- 
ment devices are always open systems that interact with external environ- 
ments. This interaction dynamically selects a preferred set of states of the 
apparatus, the so-called pointer basis. This is, in some sense, the set of 
the most stable states: if the aparatus is prepared in a pointer state, the 
interaction with the environment has a minimal effect and almost no pre- 
dictive power is lost. On the contrary, if the initial state is a superposition 
of pointer states, the interaction with the environment induces correlations 
and the state of the system tends to evolve into a mixture of pointer states. 
This process was called decoherence and this is the sense in which we will 
use this word here. Within this context, there are several important issues 
that require further attention. The most important one seems to be the defi- 
nition of an appropriate measure of stability that may be used to determine 
the pointer states (see Zurek's contribution in this conference). We will 
study a model describing a particle interacting with an environment formed 
by a collection of harmonic oscillators. In this case the pointer states of the 
particle seem to be closely related to coherent states and decoherence is the 
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process that supresses interference between coherent states (our diagnostic 
tool will be described later) . The action of the model is the following: 

t 

S[x, q]= J ds - ^qX^^ + ^ ^mniql - ^IqI) - ^nxq-n (1) 





We will assume that there are no initial correlations between the system 
and the environment (i.e. the initial density matrix factorizes) and that the 
initial state of the environment is in thermal equilibrium at temperature T. 

In Section 2 we will describe generic features of the evolution of a system 
interacting with a general environment. We will study in detail the case in 
which the initial state is a superposition of two coherent states. In Section 
3 we will illustrate the fact that decoherence strongly depends on some 
properties of the environment. We will also illustrate in what sense position 
is a preferred observable in the model (where position eigenstates are not 
pointer states). In the Appendix we outline a simple derivation of the master 
equation for a general environment. 

General Properties of the Reduced Dynamics 

Due to the interaction with the environment, the evolution of the system 
is non-unitary since it is affected by a stochastic noise and a "dissipative" 
force (the word dissipation is used here in a rather vague sense). Noise and 
dissipation, are entirely determined by two properties of the environment: 
the spectral density /(a;) and the initial temperature. The spectral density, 
defined as /(a;) = Ylm ^('^~'^n)C'n/2™n'^n) characterizes the number density 
of oscillators in the environment and the strength of their coupling with the 
system. Therefore, in order to analyze how decoherence depends on the 
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environment we can study how this process changes when varying I{u}) 
and the temperature since these are the only two environmental properties 
"seen" by the system. 

Luckily enough, the reduced dynamics has some very general features 
that are entirely independent of the spectral density and the temperature. 
One of the most striking and important ones is the fact that the reduced 
density matrix always satisfies a master equation that can be written as 
follows (we use ^ = 1): 
i dtpr{x,x',t) = < x\[Hrenit), pr]\x' > -i^{t){x - x'){dx - dx>) Pr{x,x',t) 

-iD(t){x-x'f pr{x,x',t) + f{t){x-x'){dx + dx>) Pr{x,x',t) 

(2) 

This equation depends on the spectral density and the temperature only 
through the coefficients appearing in the right hand side: the physical fre- 
quency entering in Hrenit), the friction coefficient j{t) and the diffusion 
coefficients D{t) and f{t) are time dependent functions that vanish initially 
and depend on the environment in a fairly complicated way. 

The validity of (2) for a general environment at arbitrary temperature 
has been recently demonstrated by Hu et al. (1992) and is an interesting dis- 
covery that generalizes previous results concerning the nature of the master 
equation (see Unruh and Zurek (1989), Caldeira and Leggett (1985), Haake 
and Reiboldt (1985)). The result is also surprising since a general envi- 
ronment generates a non-Markovian evolution for which one expects highly 
nonlocal integral kernels in the master equation. However, for this model 
the non-Markovian effects can be fully encoded in the time dependence of 
the coefficients. The interested reader can find a simple proof of equation 
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(2) (the simplest I could think off) in Appendix 1. 

Equation (2) is a very useful tool to study generic properties of the 
evolution and can be exactly solved for some simple initial conditions. To 
study decoherence we will consider the following initial superposition of 
coherent states: 

iif{x,t = 0) = Nexp{-^^^^+iPox) + iVexp(-^^^^ - iPQx) (3) 

where N is a constant. In this case it is possible to solve the master equation 
and show that the Wigner function constructed from the reduced density 
matrix is: 

W{x,p,t) = Wi{x,p,t) + W2{x,p,t) + W,nt{x,p,t) (4) 

where 

W, (x, p, = ^ exp _ slip ^pe- (3{x T xc)?) 



5 x^ 

Wint{x,p,t) = 2N'^ — exp(--p- - 5l{p - (3xf) cos((/)pp + {(px - f3(pp)x) exp(-A,^i) 

(5) 

The functions Xc{t)., Pc{t), Si{t), f3{t) ,(f)x{t).,(l)p{t) and Aini{t) depend on 

2 

the environment (and on the constants Lq, Pq and 6 that appear in (3)) 
in a rather complicated way. For the sake of brevity, we will not discuss 
here the behavior of all these functions (see Paz et al (1992) for details) but 
concentrate on A^^^ which is the only one relevant for decoherence. Thus, to 
quantify the importance of interference at a given time we will use the peak 
to peak ratio between the interference and the direct terms in the Wigner 
function, a quantity closely related to Ai^t- 

. , ^ 1 Wint{x,p)\peak ... 
ex.p{-Aint) = (6) 



(Wi{x,p)\peakW2{x,p)\peak) 



As the two initial wave packets have a finite overlap, the above function 
satisfies A^^^ < 5^Pq+Lq/5'^. The system decoheres when A^^t irreversibly 
grows to a value that is large with respect to unity (which can only occur if 
the initial peaks are well separated, i.e. delta^PQ + » 1). 

To analyze how the evolution of Aj^^ is affected by the environment it 
is convenient to use the master equation (2) to show the following identity: 

Air,t = D{t)4 - 2f{t)M^:c - f3<t>p) (7) 

The first term carries the effect of normal diffusion and always produces 
decoherence since increases the value of On the contrary, the sign of 

the second term in (7) may vary in time depending upon the relation between 
(f)p and (f)x- Equation (7) can be approximately solved if one neglects the 
anomalous diffusion and considers D(t) as a constant, two conditions met 
by an ohmic environment in the high tempearture regime (see next section) . 
In this case it can be shown that, for an initial state with Pq = 0? ^intif) — 
4.LlDt/{l + 4.D5H) and that the "decoherence rate" is T^^^ = 4i^o-^ - 
8Lgm7oA;5T (see Paz et al. (1992)). However, this solution is no longer valid 
when one moves away from the ohmic environment in the high temperature 
regime or when considers more general initial states. In fact, the behavior of 
A^nt strongly depends on the initial conditions (that enter into (7) through 
the functions (px and (pp whose initial data are (f)x = Pq, 4>p = Lq) and 
decoherence will be drastically different in the case Lq = 0, Pq ^ where 
the two initial gaussian are spatially separated than when Pq = and Lq ^ 
(where the coherent states are separated in momentum). 

Decoherence and the environment. 
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A wide and interesting class of environments is defined by a spectral 
density of the form I{u) = exp (— ^) where A is a high frequency 

cutoff and n is an index that characterizes different environments. We will 
consider two examples: n = 1 which is the largely studied ohmic environ- 
ment (Caldeira and Leggett (1985)) and n — 3 which is a supra-ohmic 
environment used to model the interaction between defects and phonons in 
metals, Grabert et al (1988), and also to mimic the interaction between a 
charge and its own electromagnetic field, Barone and Caldeira (1991). 

Using these two environments we want to illustrate how strongly deco- 
herence depends on the spectral density. It can be shown that the process 
is much more inefficient in the supraohmic than in the ohmic case because 
the final value of the diffusion coefficient D{t) is much smaller in the former 
than in the latter environment (as n = 3 corresponds to a bath of oscilla- 
tors with an infrared sector substantially weaker than n = 1, the dissipative 
and diffusive effects are expected to be weaker). The time dependence of 
the diffusion coefficient for these two environments has been described by 
Hu et al (1992) and has a rather generic feature: D(t) vanishes initially 
and develops a very strong peak in a time scale of the order of the col- 
lision time Tyv = 1/A. Its value after the initial peak is D{t\) ~ tti^qA, 
approximately the same for all environments. After this initial cutoff dom- 
inated regime, D{t) approaches (in a dynamical time scale) an asymptotic 
value that depends on the environment (in the high temperature regime, 
D{t) —>■ 2m^QkBT for n = 1 while vanishes for n = 3). Thus, there is no 
generic long time behavior but a quite universal short time regime. One 
may thus wonder if this general initial behavior produces a rather universal 
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decoherence. We will argue here that this is not the case. The impact of the 
initial peak has been analyzed in detail (see Unruh and Zurek (1989)) and 
it was shown that in some cases may completely wipe out the interference 
effects. However, the physical significance of the decoherence produced by 
the initial peak is rather questionable since this jolt is certainly related to 
the initial conditions that do not contain correlations between the system 
and its environment. In fact, such correlations are likely to wash out the 
initial peak, Grabert et al (1988). 

To discredit even more the role of the initial peak on decoherence we 
would like to point out that its efi'ect can be made completely inocuous by 
appropriately choosing the initial conditions for the system. This is well 
illustrated by the supra-ohmic environment where the asymptotic value of 
the diffusion coefficient is too small to produce decoherence and all the 
effect, if any, should come from the initial peak. In Figure 1 we plotted Ai^t 
for the ohmic and supraohmic environments. We considered an harmonic 
oscillator with renormalized frequency Qr and fixed 70 = O.SQr, A = bOOQr 
and k^T = 25000f2r (high temperature regime). 

We can notice that in the ohmic environment decoherence is very fast. 
For the initial condition I (Lg = 35, Pq = 0) it takes place in a time of the 
order of Tyv while for condition II {Lq = 0, Pq = 3/6) it requires a time 
that is also much smaller than On the other hand, in the supraohmic 

environment of Figure l.b decoherence goes as in the ohmic case for con- 
dition I while no net decoherence is achieved for condition II. In this case 
the initial growth of ^j^^ is followed by a plateau and a decreasing regime 
during which coherence is recovered! The reason for the drastic difference 
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between the fate of conditions I and II in the supraohmic environment is 
clear: decoherence can only be produced by the initial peak but the interac- 
tion between the system and the environment is initially effective only if the 
two coherent states are spatially separated. The non monotonic behavior of 
Ajnt seen in curve (II) of Figure l.b is due to the anomalous diffusion that 
cannot produce any net decoherence since the sign of the second term in the 
r.h.s. of (7) changes with time. 

The above example not only illustrates the strong dependence of deco- 
herence on the spectral density but also clarifies in what sense position is 
an observable that is preferred by the interaction. In fact, in general, co- 
herent states that are spatially separated decohere much faster than those 
separated only in momentum (see Paz et al (1992) for more details). 
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Appendix: Derivation of the Master Equation 

I outline here a simple derivation of the master equation (2) based 
on the properties of the evolution operator of the reduced density matrix. 
This propagator, which we denote as J(t, to) and is defined so as to satisfy 
Predi^) = ^o)Pre(i(^o); has a path integral representation of the following 
form: 



J{x,x',t I XQ,XQ,tQ) = / Dx / Dx exp — I S[x] — S[x] > F[x,x] (8) 



where F(x,x') is the Feynman- Vernon influence functional that arises due 



X 



x' 
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to the integration of the environment variables. For the model we are con- 
sidering, this functional is well known and can be written as (see Grabert 
et al (1988)): 

t s 

i ln{F[x,y]) = J ds J ds'{x-y){s) ri{s - s'){x + y){s') - ii^{s - s'){x - y){s')^ 


where vis) and r]{s) are the noise and dissipation kernels defined in terms 
of the spectral density: 

(jj 



r)(s) = — / du}I(u>) sm(u)s), vis)— / dwliuj) coth.(— — — )cos(a;s) 
Jo Jo ^i^bT 

As the integrand of (A.l) is gaussian, the integral can be exactly computed 
and the result is (written in terms of the variables ^ = x — x' ^ X = x + x'): 

J{X, t; Xq, Co, ^o) = ^ exp{ibiX^ + 262^0^ - ^^'S^^O - ibAXo^o) X ^ ^ 

^TT (9) 

X exp(-aii^2 - au^^o - 022^0) 
where the functions b}^{t) and a^j(t) depend on the environment and can be 
constructed in terms of solutions to the equation: 

rs 

u{s) + nlu{s) + 2 ds'r]{s - s') u{s') = (10) 

Jo 

Thus, if ui and U2 are two solutions of (A. 3) that satisfy the boundary 
conditions ui(0) = U2{t) — 1 and ui{t) — ^2(0) = we can write: 



2 bi{t) =U2{t), 2 63(t) = {^2(0), 2 h2{t) = ui{t), 2 h^it) = iiiiQ) 

»i rt (11) 



= 1 \ 1 I ds ds'ui{s) Uj{s') v{s - s) 
1 + Oil Jo Jo 



The derivation of the master equation can be done as follows by simply 
using equations (A. 2) and (A. 4): Let us take the time derivative of (A. 2) 
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and write 



.... 2 2 

(12) 

If we multiply (A. 5) by an arbitrary initial density matrix and integrate over 
the initial coordinates ^0 cind Xq, we will obtain an equation whose left hand 
side is pr{x, x', t). In the right hand side we will find terms proportional to 
p{x, x', t) that look like some of the ones appearing in the right hand side of 
equation (2). The only potentially problematic terms are the ones that in 
(A. 5) are proportional to the initial coordinates ^0 Xq. However, their 
contribution can be easily shown to be local by realizing that the propagator 
J(t, to) satisfies: 



eo J(x, e, t; Xq, Co, to) =(^e + f %) J(x, e, t-, Xo, Co, to) 
XoJ{x, t; xo, Co, ^o) H-x^r - r^i - + tt^)^ + rf^^)'^^^' ^' ^o) 

02 02 02 0203 0302 



Using these equations, the right hand side of (A. 5) can be written in terms 
of the final coordinates X and ^ and the derivatives with respect to them, 
which implies that the master equation is local. To show that this equation 
is given by (2), we just have to demonstrate that the coefficients associated 
to terms like 6^ or Xdx cancel and this can be done by exploiting some 
general properties of the coefficients hj. and a^j that follow directly from 
their definition in (A. 4). In fact, using relations such as 022 = —^4012/^2 
(whose proof we omit), it is possible to show that the coefficients of equation 
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(2) are: 




h); 



lit) 




ren 



an - 4aii6i + di2— - T-(2aii + a^—) 
ai2 _ ^2^12 _ ^^^^ 
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Figure Captions 

Figure 1: The degree of decoherence A^^t is plotted as a function of time 
(which is measured in units of ^r^) for the ohmic (a) and supra-ohmic (b) 
environments in the high temperature regime. Curve (I) corresponds to 
an initial condition where the initial coherent states are spatially separated 
(Lq = 3S, Pq = 0) while for curve (II) the initial separation is in momentum 
(Lo = 0, Po = 3/5). 
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Questions 

Unruh: Your master equation is local in time. For an arbitrary spectral 
density I would strongly expect that the equations are strongly non local in 
time. Why aren't yours? 

Paz: There are two observations one can intuitively make for the model 
described by equation (1). On the one hand we expect that a general envi- 
ronment will produce non-Markovian effects and that the master equation 
will be non local in time. On the other hand, the (reduced) evolution oper- 
ator must be gaussian since the problem is linear. The crucial observation 
is that if one admits a gaussian evolution operator, the master equation is 
always local provided the matrix mixing "old" and "new" coordinates in the 
propagator can be inverted. Taking this into account, it is surprising to me 
that the existence of a local master equation has not been noticed until so 
recently. 

Morikawa: Why do you get a local coefficient 7(t)? 

Paz: I refer to the answer I gave to Prof. Unruh's question. 
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